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Using the Green's function formahsm, an ab initio theory for band structures of crystals is de- 
rived starting from the Hartree-Fock approximation. It is based on the algebraic diagrammatic 
construction scheme for the self-energy which is formulated for crystal orbitals (CO- ADC). In this 
approach, the poles of the Green's function are determined by solving a suitable Hermitian eigen- 
value problem. The method is not only applicable to the outer valence and conduction bands, it 
is also stable for inner valence bands where strong electron correlations are effective. The key to 
the proposed scheme is to evaluate the self-energy in terms of Wannier orbitals before transforming 
it to a crystal momentum representation. Exploiting the fact that electron correlations are mainly 
local, one can truncate the lattice summations by an appropriate configuration selection scheme. 
This yields a flat configuration space; i.e., its size scales only linearly with the number of atoms 
per unit cell for large systems and, under certain conditions, the computational effort to determine 
band structures also scales linearly. As a first application of the new formalism, a lithium fluoride 
crystal has been chosen. A minimal basis set description is studied, and a satisfactory agreement 
with previous theoretical and experimental results for the fundamental band gap and the width of 
the F 2p valence band complex is obtained. 

PACS numbers: 71.10.-w, 71.15.Qe, 71.20.-b 



I. INTRODUCTION 

Band structures are of fundamental interest to the 
solid-state physicist as they reveal important properties 
of crystalsiiiSiiii^ The determination of band structures 
on ah initio level is therefore an important issue in the- 
oretical solid state physics. Nowadays predominantly 
density functional theory (DFT)£i& with the local den- 
sity approximation (LDA) or its variants arc employed to 
treat this problem. ^'"^ DFT focuses on the ground-state 
electron density and yields ground-state properties. In 
several cases, the interpretation of the Kohn-Sham or- 
bital energies as quasiparticle energies also turns out to 
be successful and particular attention has been paid to 
the LDA, due to its numerical simplicity. However, in in- 
sulating materials, the LDA tends to underestimate the 
band gap, e.g., Ref.|3. 

As is well known, it is difficult to improve band struc- 
ture calculations based on density functional theory. This 
is in particular the case when electron correlations are 
strong such that in quite some casesi*i2*ii*i2iiii4 the one- 
electron picture becomes inappropriate and corrections 
are needed. Nevertheless, a number of ingenious meth- 
ods have been devised to improve band structure results 
based on density functional theory. For example, the cal- 
culations have been supplemented by a GW treatment, 
e.g., Refs Bfini. or the so-called LDA-hC/ methodii 
Other improvements concern the use of optimized effec- 
tive potentials'^ such as the exact exchange potentiali^ 
instead of its local, simplified form. Time-dependent 
DFT is another route to treat excited states more rigor- 
ously. Furthermore, LDA calculations have been coupled 
with various extensions of the coherent potential approx- 
imation (CPA) like the dynamical CPASSiSi or dynamical 
mean field theory (DMFT),^'^i^3,24,25,26,27 ^isp Faddeev's 



method^- of treating the three-particle i-matrix has been 
applied in conjunction with LDA band structure calcu- 
lationsi^S The same holds true for different forms of the 
projection operator techniquei^ Yet one must face the 
fact that the various approximations remain uncontrolled 
and therefore have to be reconsidered from case to case. 

When one is aiming at controlled approximations, two 
different routes offer themselves. One is bound to quan- 
tum Monte Carlo (QMC) calculations, e.g., Ref. Isol 
This is a simple and straightforward way of dealing with 
the many-body problem. Yet it suffers from one seri- 
ous shortcoming, the so-called sign problem when treat- 
ing fermions. Therefore, we are pursuing the second 
route, ft starts from Hartree-Fock band structures and 
improves them by including correlations by means of 
many-particle techniques similar to the ones employed 
for molecules.''Q-^'-^'--^-'^-'^'-'^5i36,37,38,39,4Q,4i Thereby one 

can take advantage of program packages like wannier 
which has been written by Shukla et alMi^ or crystal 
which originates from Torino.'^'-- Both provide Hartree- 
Fock bands as well as Wannier orbitals for crystals. These 
orbitals are especially suited for an ah initio treatment of 
electron correlations because of their local character 

A number of correlation calculations have 
been performedi^iSiSi^ A local Hamiltonian ap- 
proach22i22i2i^Si22i^ was shown to improve Hartree-Fock 
energy bands substantially. Suhai used Toyozawa's elec- 
tronic polaron model^ to repartition MP2 pair energies 
to estimate quasiparticle band structures i2iiS4SiiS4L^ 
By inserting the orbital energies into the energy- 
dependent self-energy, this model was shown to be 
a special case of the outer valence Green's func- 
tions (OVGFs)i2i^2*^ that were derived in terms of 
crystal orbitals by Liegeneri2iiS42iS By applying 
the approximation to the self-energy of Igarashi et 
al.^S^^ quasiparticle band structures were obtained by 
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Albrecht et aZ.ffl-S 

In the present paper, we derive the algebraic diagram- 
matic construction (ADC) scheme^^iS^i^ for crystals. 
The ADC scheme has proven to be superior to the OVGF 
method in molecular studies)^ and numerous works 
have been carried out over the last two decades, e.g., 
Refs.[T^13''l4'. including studies of oligomers and clusters 
chosen to model infinite chains or crystals. ^^■^^"'^^ The 
ADC scheme is a method to approximate the Feynman- 
Dyson perturbation series for the self-energy and contains 
sums of certain proper and improper diagrams to infinite 
orderiSiSSiSi The basic properties of the ADC scheme 
and their derivations, as applied to molecules, are found 
in Rcfs. 62 63 64. Among those sizc-cxtensivity is an im- 
portant one since it is crucial when solids are considered. 
Furthermore, the ADC method is known to be robust and 
facilitates also to study strong electron correlations due 
to the efficient and stable evaluation of the one-particle 
Green's function in terms of a Hermitian eigenvalue prob- 
lem. In molecules strong correlations for instance occur 
when inner valence electrons are treated >i2iiii4 With the 
crystal orbital ADC (CO-ADC) devised here, the spec- 
tral representation of the one-particle Green's function is 
obtained. We suggest a formulation of the expressions 
in Wannier representation instead of crystal momentum 
representation in order to be able to exploit the local 
character of the correlation hole around an electron es- 
pecially in the case of crystals with a large unit cell. 

As a first application of the CO-ADC method to 
a three-dimensional crystalline solid, we have chosen 
lithium fluoride which occurs in nature as the mineral 
griceite. LiF not only has a wide range of technologi- 
cal applications like in x-ray monochromators or in fil- 
ters for ultraviolet radiation, e.g., Ref. |6^ (and Refer- 
ences, therein), but is also interesting for a number 
of fundamental physical reasons. It is considered to be 
the "most ionic substance" and a prototypical insulator 
which manifests in its very large fundamental band gap 
of 2 ± 0.1 eV that is the largest one found in na- 

ture apart from exotic systems. Some authors even con- 
sider LiF to be comparable to a He-Ne rare-gas solidiS^ 
Its optical spectra show strong excitonic effects, which 
complicates the experimental determination of the band 
gap^^ by optical spectroscopy. Similarly, many-particle 
effects have to be accounted for in the measurement of 
the widths of the F 2p and F 2s valence bands. 

LiF has been thoroughly studied both experimentally 

and theoretically|38,39i43,65,67,68,69,70,71,72,73 pooig 

survey early experimental and theoretical dataiSMS Re- 
cent studies of LiF comprise density functional theory 
calculations at the LDA level supplemented by an in- 
clusion of correlation effects in terms of the GW approxi- 
mationJ^'®'^''''^ Ab initio investigations comprise Hartree- 
Fock studies of ground-state properties^SiZi and of the 
band struct urei^24^2ii& An accurate treatment of electron 
correlations for the ground-state properties of LiF (and 
other alkali halides) has been carried out by Doll and 
StoU.^^ Quasiparticle band structures have been obtained 



by Kunz''° on the basis of the electronic polaron model^ 
and by Albrecht^ with the Green's function approach of 
Igarashi et al.m^ 

The article is structured as follows: SectionllTlputs for- 
ward an ab initio foundation for band structures based on 
the one-particle Green's function. In Secs. lIIII and llVI we 
derive the CO- ADC approximation to the self-energy in 
terms of Wannier orbitals. SectionlVlcombines the results 
of the previous sections to formulate the band structure 
problem as a Hermitian matrix eigenvalue problem. Cri- 
teria for the truncation of lattice sums are discussed in 
Sec. EJ In Sec. |VlII the apphcation of CO- ADC theory 
to lithium fluoride is described. Conclusions are drawn 
in Sec. lVml 



II. BAND STRUCTURES 



The nonrelativistic Hamiltonian of a crystal in atomic 
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where flxed nuclei are assumed. Here N is the number 
of electrons in the crystal, Nq denotes the number of 
unit cells, and M is the number of nuclei per unit cell. 



RiA 



Za stands for the charge of nucleus A in unit 
cell Ri, and {t^u ^ rj^ . is the distance between the n-th 

electron and the A-th nucleus in unit cell Rj. Finally, 
\frn — 'Tnl rcprcscnts the distance between the m-th and 
n-th electron and I^jI j^ ~ ^j^.gl denotes the distance be- 
tween the nuclei A, B of charge Zj^^ = Z^, Zj^^ = Zb 
in unit cells Rj. 

The Schrodinger equation with the Hamiltonian 
can be solved for the ground state in a restricted, closed- 
shell Hartree-Fock approximationfifi*ii*£^iS*£& yielding a 
set of spin Bloch orbitals where k is the crystal 

momentum and p is a band indexi-^i^i'^i'^i^i-^°i-^-^i'''^i'^^i'^^ We 
use these orbitals to represent the Hamiltonian ^ in 
second quantization, yielding the so-called Bloch or crys- 
tal momentum representation."^ It is decomposed accord- 
ing to M0ller and Plesset2i22i2S into the Hartree-Fock 
part -ff^F ^n'i the residual interaction H^^, a suitable 
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form to apply perturbation theory: 
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The energy bands (or Bloch orbital energies) are denoted 
by e^p, and the operators ct create (annihilate) 

electrons in ipftpir). The negative of the Hartree-Fock 
potential is given by W^g^ ^ ^ = - E l^fc, p fe. [£ 
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the occupation numbers n^^ being unity for occupied 
and zero for unoccupied (or virtual) Bloch orbitals (and 
"fcp ~ -'^ ~ "fcp)- Generally, we use the band indices i, 
j , m, n, . . . to denote occupied Hartree-Fock bands, a, 
6, c, d, ... for unoccupied bands, and p, q, r, s, . . . 
for bands which are occupied or unoccupied. The two- 
electron integrals in Eq. Q are defined with respect to 
Bloch orbitals byT7i7?i7? 
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The antisymmetrized two-electron integrals are 
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and consist of a Coulomb term Vr , r „ r ^ r . 



(4) 
and an 

exchange term V^^ k^sk r- crystal momentum 
vectors in Eqs. 12}, ©, and (@J) are restricted to the first 
Brillouin zone. Furthermore, we assume that the first 
Brillouin zone is discretized by the Born von Karman 
boundary conditions, or equivalently, that numerical in- 
tegrations in reciprocal space are carried out in terms of 
a Monkhorst-Pack net.^'^'^" In both cases, integrations 
over the Brillouin zone are approximated by finite sums. 
The one-particle Green's function (or particle propa- 



gator) in terms of Bloch orbitals 
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with Wick's time-ordering operator T. \ ) is the exact 
ground state of the N particle system. Due to transla- 
tional symmetry of the Hamiltonian the one-particle 
Green's function depends only on one crystal momen- 
tum^i.e., G^^^, ^{t,t') ^ 5^^^, Gpq{k,t,t'). 

Similarly, the time independence of the Hamilto- 
nian (Q) causes the particle propagator © to depend 
only on the time difference t — ^'.giigSiS Fourier trans- 
forming Gpq{k,t — t') with respect to t — t' yields the 



one-particle Green's function in energy space Gpq{k,uj) 
which can be recast in terms of the spectral or Lehmann 
representation asSLSSiS 
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Here n numerates the excited states of the ± 1 par- 
ticle system. The negative of the pole positions in 
Eq. (|6a|l is given by either the electron affinities An{k) = 
Eq — E^^^{k) or the ionization potentials /„(fc) — 
E^-^{-k) - E^ where £;^=^i(±fc) are the energies of 
the different excited states of the ± 1 particle system. 
The summands ztir; are required to ensure the conver- 
gence of the Fourier transformation. Whenever such a 
factor occurs, we implicitly take the limit 7y ^ 0+ iSLSiS 
The pole strengths in Eq. are given by the transi- 
tion amplitudesSi 



k p 



^^+\k)) , 



(fc) = (*rH-fc)l%pl*o^> 
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where |^'^^^(±fc)) denote excited states of the N ± 
1 particle system with crystal momenta ±fc. The pole 
strengths can be interpreted in terms of the spectral 
intensities observed in photoelectron spectroscopy ex- 
periments, similarly to the molecular case discussed in 
Ref.H^ 

To obtain electron affinities and ionization potentials 
from the band structure of three-dimensional crystals, 
one has to specify the energy of the added (removed) 
electron at the Fermi level of the neutral crystal. ^^'^^ This 
is usually done by introducing chemical potentials /i^ 
which are added to the pole positions of the one-particle 
Green's function, i.e., A„(fc) -I- and /„(fc) -I- ^~ . 

The (fundamental) band gap is the smallest difference 
between the energies for removing an electron from and 
attaching an electron to an N particle system: 



E, 



= I - - 



(8) 



Here A stands for the largest electron affinity and / des- 
ignates its smallest ionization potential. Correspond- 
ingly, E^^^ denote the energies of the electron attach- 
ment (removal) state with the lowest (highest) energy 
and Eq represents the ground-state energy. 

The self-energy S(fc,w) with respect to the resid- 
ual interaction is defined by the Dyson equa- 

^^.-,^ ^3.55.81.82.83 

G(fc, oj) = G°{k, uj) + G°{k, oj) S(fc, uj) G{k, uj) (9a) 
= G"(fc, uj) + G"{k, uj) S(fc, uj) G°{k, uj) (9b) 
... , 



4 



which can be solved formally by 

G{k,uj) = [G°{k,uj)-^ - S(fc,cj)]-i 

where 



k p, k q 
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is the free Green's function with respect to the M0ller- 
Plesset partition jSJ- In fact, for crystals with a band 
gap, the occupation numbers in Eq. are independent 
of k, i.e., — Up and = Hp. 

The self-energy can be decomposed into an uj indepen- 
dent part, the static self-energy I]°°(fc), and an ui depen- 
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dent part, the dynamic self-energy M{k,uj'^ 50,55,87 



S(fc,w) = i:°°{k) + M{k,uj) 



(12) 



where lim M{k,Lu) = 0. The dynamic self-energy is 

UJ — >±oo ' 

considered first in the ensuing Sec. IIIII as the static self- 
energy is determined using the former one in Sec. IIVI 



III. DYNAMIC SELF-ENERGY 

A. Crystal momentum representation 

The partition |5J of the Hamiltonian has the same form 
and meaning as the M0ller-Plesset partition for canoni- 
cal molecular orbitalsi^Ii'^S''''^, used frequently in molecu- 
lar physics. Therefore, all equations derived on the ba- 
sis of the molecular M0ller-Plesset partition are in com- 
plete analogy to the equations in the crystalline case. 
The only difference between molecules and crystals is the 
occurrence of composite indices in the Hamiltonian |(2Jl, 
consisting of a crystal momentum vector k and a band 
index p which vary independently. Given an equation 
in terms of molecular orbitals, it can be written immedi- 
ately in terms of Bloch orbitals by replacing all molecular 
orbital indices by composite indices. Afterwards one can 
exploit (translational) symmetry to simplify the equa- 
tion. 

The dynamic self-energy M{k,Lu) is represented in 
terms of the 2pl/i/2ft,lp-propagator (two-particle-one- 
hole/two- hole-one-particle-propagator) and possesses the 
spectral representation^^iSLS 
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where (^ni^^) denote the pole positions of the 2plh/2hlp- 
propagator and the m^'^^\k) are termed Dyson ampli- 
tudes. The retarded dynamic self-energy M^^{k,uj) and 



the advanced dynamic self-energy Mp^{k,uj) are associ- 
ated with excitations of the ± 1 particle system, re- 
spectively. Yet the pole positions 0Jn{k) do not directly 
correspond to the energies of physical states as the (dy- 
namic) self-energy is only defined in conjunction with the 
Dyson equation 

The algebraic diagrammatic construction scheme is a 
stable and efhcient method to determine the spectral rep- 
resentation of the dynamic self-energy H13() . The con- 
struction starts by making the ansatz 

Af±(fc,c.) = C/±t(fc) (,, 1 _ K^{k) - C^{k))-^ U^{k) , 

(14) 

which is termed general algebraic form or ADC form. 
The vectors U^{k) are called modified coupling ampli- 
tudes for the crystal orbital kq. Both K^{k) and C^{k) 
are Hermitian matrices, C^{k) being referred to as mod- 
ified interaction matrices^ while K^(k) are assumed to 
be diagonal. Explicit expressions for these matrices are 
given below in Eq. p8(l . 

The ADC form (|14|) reproduces the analytic structure 
of the Feynman-Dyson perturbation series for the dy- 
namic self-energy as can be seen from its expansion into 
a geometric series: 
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(15) 



The n-th order approximation to the ADC form 114|l is 
constructed by inserting the perturbation expansions 

U^{k) = U^'^^\k) + U^'-^\k) + ... , 
C±(fc) = C^^^\k) + C^^^\k) + ... , 



(16) 



in terms of the residual interaction iJ^f into Eq. I|15|l . 
and comparing the resulting expressions with the dia- 
grammatic expansion of the dynamic self-energy up to 
n-th order. This constitutes a scheme which is denoted 
by ADC(n). Note that the analytic structure imposed by 
Eq. (|15|) occasionally requires one to associate linear com- 
binations of the analytic expression of several diagrams 
with particular terms in the above expansion. 

In ADC(2) and ADC(3), each entry of the C7±(fc) and 

the K"^ (k) + (k) in Eq. ((Tl|l is characterized by one or 
two, respectively, arrangements of two particles and one 
hole (2plh) for M^q{k,uj) and two holes and one parti- 
cle {2hlp) for M~^{k,uj). In brief notation these arrange- 
ments are referred to as 2p\h- and 2/ilp-configurations. 
All such configurations that can be formed with the 
one-particle basis set underlying a Hartree-Fock calcu- 
lation constitute the configuration space in ADC(2) and 
ADC(3). The configuration space is enlarged in ADC(4) 
by 'ip2h- and 3/i2p-configurationsi^In general, every sec- 
ond order in ADC, the configuration space enlarges by 
the next higher excitation class. Restricting our discus- 
sion to 2plh- and 2/ilp-configurations, we can rewrite 
Eq. ini) as 
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The prime on the summation symbol indicates that the 
summation only runs over indices with lia < /2 6 and 
I'l a' < I2 h' to avoid double counting of contributions. 
The corresponding ansatz for M~{k,u>) is formally iden- 
tical to Eq. H17() apart from the changed occupation num- 
bers of the band indices i —i- a, a i, b j, i' a', 
a' — > i' , and b' — > j'. The summation variables h, I2, 
I'l, and ^2 are crystal momentum vectors like k. The 



quantities M^{k,Lu), U^{k), K^{k), and C^{k) are in- 
variant under translations by an arbitrary lattice vector, 
a property which is taken into account for the internal 
summations in the ADC form p7|l . 

Explicit molecular ADC equations up to fourth or- 
der are given in Ref. Here we show only the ADC 
equations for the dynamic self-energy in terms of crystal 
(Bloch) orbitals up to second order [CO- ADC (2)]: 
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(18a) 
(18b) 

(18c) 



The equations for U (k), K (k), and C (k) are for- 
mally identical to Eq. 118|) apart from the change of la- 
bels as for the adaptation of Eq. H17|l to MzSk^uj) and 



the occupation numbers in Eqs. H18a|l and Ijl8b(l which 
are in this case n;- _|_;- ".f^ j ".f^ j ■ Expression (|18|l is de- 
rived in Ref. |8^ utilizing a Gaussian basis set expansion 
of the Bloch orbitals. 

The spectral representation (|13|l of the dynamic self- 
energy is obtained from the ADC form (|14|l by solving 
the Hermitian eigenvalue problem 



{K^{k) + C^{k)) Y^{k) = Y^{k) n^{k) 



Y^^k)Y^{k) = 1 



(19) 



The vector notation in Eqs. (|14|l and H19|) is a convenient 
means to sum over all intermediate crystal momenta and 
band indices. The diagonal matrix ft^{k) contains the 
eigenvalues of the secular matrix K^{k) + C^{k), and 
Y^{k) denotes the eigenvector matrix. The eigenvalues 
are the pole positions of the dynamic self-energy, i.e., 
(^) — ^nn(^); whilc thc Dysou amplitudes in Eq. 113(1 
are obtained via 



m±'(")(fc) = [7p±t(fc)y±.(«)(fc) 



(20) 



where the n-th column vector Y^'^'^''{k) of Y^{k) is used 
here together with the adjoint of U^{k). Note that. 



unlike CO- ADC (2), in CO- ADC (3) (and higher orders) 
the eigenvectors y^'("'(fc) couple several 2plh-, 2hlp- or 
higher excited configurations via the modified coupling 
amplitudes H18c|l which are nonzero in this case. 



B. Wannier representation 

We derive in the following an expansion of the dy- 
namic self-energy in terms of generalized Wannier or- 
bitals; such orbitals are localized within unit cells and 
allow one to apply cutoff criteria inside a cell and be- 
tween a cell and other cells. This affords one the ability 
to apply fine-grained configuration selection which is in- 
dependent of the actual choice of the unit cell. In partic- 
ular, it enables the treatment of crystals with large unit 
cells, not amenable to Bloch-orbital-based approaches 
[Sec.|VlliWi,35 

The generalized Wannier transformationii2i2i*2^i2242ii2^ 

and its inverse which mediate between Bloch and Wan- 
nier orbitals are 
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^^Z^,,(fc)e-''=«^,,^(f),(21a) 
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where K denotes the number of bands involved in the 
transformation and I4{k) is a suitable unitary matrix 
as generated, e.g., by the locahzation procedures in 
Refs. I9li92l. Both Bloch and Wannier orbitals are nor- 
malized to unity upon integration over the entire crystal 
which consists of iVo unit cells. Generally, we denote the 
Wannier orbital indices that result from occupied bands 
with K, X, fj,, ly, . . . , while the Wannier orbital indices of 
unoccupied bands are denoted by a, 7, (5, . . . . Indices 
of Wannier orbitals from occupied or unoccupied bands 
are referred to by p, a, t, v, . . . . 

For an evaluation of the one-particle Green's function 
in terms of the Feynman-Dyson perturbation series, it is 
essential that there is no mixing between occupied and 
virtual Wannier orbitals; otherwise Wick's theorem can 
no longer be appliediSLSSiS^ To fulfill this requirement, 
the Wannier transformation (|21a|) is applied to the occu- 
pied and virtual Bloch orbitals separately, yielding two 
independent unitary matrices, one for the occupied bands 
and one for the virtual bands, respectively. Hence, U{k) 
in Eq. (|21|l is block-diagonal. 

On the one hand, we want to evaluate the CO-ADC 
equations entirely in terms of Wannier orbitals, in the 
so-called Wannier or local representation^^ on the other 
hand, band structures are defined with respect to the 
crystal momentum quantum number. In what follows 
three different transformation schemes to switch between 
the local representation and the crystal momentum rep- 
resentation are applied to the CO-ADC formalism and 
their physical and methodological implications are dis- 
cussed. 

The dynamic self-energy Mj^^(fc,w) = -. (lu) de- 

pends on two external Bloch orbitals ipt (r) and tp^^(r). 
Hence, carrying out the inverse Wannier transforma- 
tion l|21bl) for the external orbitals in Eq. 117|l leads to 



RM' 



(22) 



for the left-hand side of Eq. H17|l and a similar expression 
for the right-hand side where the internal summations li, 
I2, I'l, and P2 still run over the Brillouin zone. 

In order to obtain a representation of the dynamic self- 
energy M^j^{k,u!) entirely in terms of Wannier orbitals, 
it is important to note that, although matrix elements 
such as C/J" - - - - - , in Eq. itTTll only depend on 

three independent crystal momenta, they actually de- 
scribe quantities which depend on four Bloch orbitals, 
namely, V'fcp(^), i^l,+^2~k^{r), i'f^Jr), and ^(r) and 
thus a fourfold inverse Wannier transformation Ij21b(l has 
to be applied. This transformation of Eq. (|17|l to the 
Wannier representation yields 



Re R' cr^ ' 



y u+* 

91 K, 92 a, 93 /3, 
Si' 92 93 0' 



X 1 - K"^ - C^)-} . - „ (23) 

^ ' 91 n 92 a 9-i P ' 

g{ k' 92 a' 33 /3' 

xU± , , , 

R' cr ; 9i k' 92 " S3 P' 

exploiting the unitarity oiU{k) and the basic relationsi*^ 



R 
k 



(24a) 
(24b) 



The check accent on the quantities f/^, , and 
in Eq. H23(l indicates that the sums run over lattice vec- 
tors and Wannier orbital indices rather than crystal mo- 
menta and band indices as in Eq. (|17|l . The prime on 
the summation symbol indicates that g2Ct < gs (3 and 
92 ^' < 53 must hold to avoid double counting of con- 
tributions. The corresponding equation for Mz -. (uj) 

is formally identical to Eq. (|23|l apart from the changed 
occupation numbers k a, a k, /? — > A, k' ^ a', 
a' k', and /?' ^ A'. 

One obviously Hermitian form for the dynamic self- 
energy emerging from Eqs. 1)22(1 and ^2'A\ is given by 



M±(fc,c^) 



(25a) 



(k) , 



E 



± 

Ra 



(25b) 



The vector notation in Eq. H25|l combines the six inter- 
nal summations over the intermediate lattice vectors and 
Wannier orbital indices 51 k, 92 cx, 93 5i' k', g2cx', and 
in Eq. (|^. We refer to Eq. ((23) as the supercell 
form of CO- ADC. 

The matrix K + C in Eq. H25al) does not explic- 
itly dependend on k which implies that its eigenvalues 
comprise all pole positions, i.e., pole positions for each 
/c-point. Consequently, the resulting eigenvalues also do 

not explicitly depend on k. Multiplying with (k) 
from the left projects on the desired set of eigenvectors. 
Hence, the translational symmetry among the A^o unit 
cells of the crystal is not exploited to reduce the size of 
the configuration space which is unfavorable in conjunc- 
tion with the cutoff criteria discussed in Sec. lVIl In prac- 
tice, Eq. H25|l is applied to a molecular cluster formed by 
Nq unit cells of a crystal for which Born von Karman 
boundary conditions are enforced. As a consequence, 
only A^o discrete crystal momenta are in the Brillouin 
zone. The above given transformation from the Wannier 
representation to the crystal momentum representation 
is simple and robust and can be used in conjunction with 
almost any electron correlation method, but it suffers 
from finite-size effects due to the, usually, small number 



of unit cells which are considered in practice; see, e.g., 
Ref.lM 

The translational symmetry of the dynamic self- 
energy (|23|l can be exploited by applying a transla- 
tion by —R which removes the explicit dependence 
of Mj^^ _R' external lattice vectors such that 

it only depends on the difference of the lattice vectors 
according to 

Mnaii'M = MsaR'-R<ri^) ^ - R,^) ■ (26) 

This facilitates to remove the lattice summation over R 
and the prefactor in Eq. H22I) which results in (drop- 
ping the prime on R') 

R 

To continue, we give up the idea of an ADC form 
for Mjj^^^(a;), where the modified coupling amplitudes 
on the left and on the right are related by Hermitian 
conjugation, as is the case in Eq. 1)25(1 . and arrive at 

Q,tT 

- ± -± ^ 

x{uj1-K C )-^U„{k) 

R 

where a slightly changed definition of the transformed 

modified coupling amplitude (k) = ^/NqU ^ (k) is em- 
ployed which is indicated by a tilde accent on the quan- 
tities U = U , K = K , and C = C . Because of 

I 



its asymmetric nature, we denote Eq. (|28|) as the semi- 
transformed form. 

By exploiting the translational symmetry (|26|) . we 
remove the summation over all translationally equiv- 
alent octuples consisting of the two external Wan- 
nier orbitals (r) and w^i^{t) of the dynamic 
self-energy M ^ j^, ^{uj) and two 2plh- or 2hlp- 

configurations. We arrive at a summation over R' — R, 
g, and a in conjunction with the remaining six internal 
lattice vectors and Wannier orbital indices. This formu- 
lation thus provides a far better starting point for reduc- 
ing the number of configurations by means of the cutoff 
criteria of Sec. IVII than the supercell form. Using the 
unsymmetric ADC form (|28|l implies that the pole posi- 
tions still do not dependent explicitly on k. Moreover, 
the modified coupling amplitudes l(28b|l are, in contrast 
to Eq. I)25b|l . no longer related by Hermitian conjuga- 
tion; yet, for sure, the dynamic self-energy itself stays 
Hermitian. 

This way of exploiting translational symmetry in con- 
junction with the Wannier orbitals has already been ap- 
plied successfully before to devise ab initio electron corre- 
lation methods for crystals like the local Hamiltonian ap- 
proach2Si2i2i2Silii or the Green's function based method 
of Albrecht et al.,^ILmi 

The internal structure of the dynamic self-energy is 
not regarded in the two formulas ((25(1 and ((28(1 . Expres- 
sion (PH|l results from Eqs. (|^ and if^ by exploiting 
the translational symmetry of M^{k, uj) in Eq. 1(26(1 and 
utilizing its exclusive dependence on two external Wan- 
nier orbitals. However, the translational symmetry of the 

matrix K +C in Eq. 1(28(1 can also be exploited addi- 
tionally which already was harnessed to derive Eq. 1(17(1 
from its original molecular orbital formulation. We ob- 
tain 



M+(fc,^) = Y.^,,{k)u;Ak) E E ui;g^^^^^,{k) 



u: 



{K^ {k) + {k)) gig^Ka/^; 

Si 92 «' 0' 



Si, Si K',a',l3' 



x(a;l-K+(fc)-C'+(fc))-^^^_^^^ 

SI Si «' p' 

^ ; Si Si «' «' P' ' 



- , ,.jk) = ye^^ut . , 

; 3i 92 k' a' I3'\ > Ra - Ok' Si 



a' 13' ' 



Ee'''^ (-f^^ + C^) Ok 31 Q 32/3; 
a Rk' gl+Ra' gi+RP' 



(29a) 

(29b) 
(29c) 



where the identities = C/^, = , and = hold. The corresponding relations for M {k,uj) are 
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formally identical to Eq. H29|l apart from the changed 
occupation nmnbers k — > a, a — > k, /? — > A, k' — > a', 
cJ — > k', and — s- A'. Alternatively, we can derive 
Eq. H29|l by inserting the inverse Wannier transforma- 
tion (|21b|l into the CO-ADC equations for Bloch or- 
bitals H17(l which already are fully adapted to transla- 
tional symmetry. We refer to Eq. H29|l as the ]ully trans- 
lational symmetry adapted form. 

The modified coupling amplitudes (|29b|l are con- 
structed by considering only Wannier orbitals relative to 
an origin cell (or reference cell) where the hole with in- 
dex K is assumed to reside. The external orbital index a 
of a modified coupling amplitude can be viewed to repre- 
sent a Bloch orbital tp'^^ir)- This orbital interacts with 
the 2pl/i-configurations that are pinned with one lattice 
vector to the origin cell and extend with the two remain- 
ing lattice vectors gi and 52 over up to two different unit 
cells. Alternatively, one can use the translational sym- 
metry of the i/t - , ^, to arrive at 

R(T \ Q k' g-^ a' p' 



-ikR 



m ^ - - 

Oo- ; Rk.' gl+Ra' g^ + RP' 



(30) 

Now the external orbital is WQ^{r) which is independent 

of R and consequently Eq. is interpreted to describe 
the interaction of a Wannier orbital in the origin cell 
with 2pl/i-configurations centered around all lattice vec- 
tors R which are combined to give an intermediate 2plh- 
configuration with a total crystal momentum —k. 

The matrix + in Eq. (|29c() describes the cou- 
pling among the 2pl/i-configurations. One of the lattice 
vectors of the first triple of Wannier orbitals is pinned 
to the origin cell, and the two other lattice vectors gi 
and g2 are offsets to it. The remaining three lattice vec- 
tors belong to a Wannier orbital in an arbitrary cell R 
and two further Wannier orbitals in the cells g[ + R 
and 52 + R relative to the former. Obviously the en- 
tire Eq. (|29c() can be interpreted to describe the interac- 
tion of a 2pl/i-configuration of Wannier orbitals centered 
around the origin cell with another 2pl/i-configuration 
with crystal momentum k. Due to the full exploitation 
of translational symmetry, also of the intermediate 2p\h- 
configurations, the fully translational symmetry adapted 
form yields the smallest configuration space in conjunc- 
tion with the cutoff criteria of Sec. IVII 

In order to give explicit CO-ADC expressions for the 



matrices showing up in Eq. H29() , we transform the Hamil- 
tonian ||2Jl with the help of formula Ij21b|) to the Wannier 
representation. A similar decomposition as in Eq. ((21 is 
chosen for a subsequent application of perturbation the- 
ory: 



H 

ttWF 
^0 



WF _ 



V''^Re^''Re^Re ' 

Rg 

^ S_ ^Ri gR2cr g ^R2 <T 
Hi £>, R2 (y 

+ 5. ^ 



(31) 



Rl Q, R2 tr, 



Rl Q R.2<T Ri"" R3T ^ 



with eg^ 



F- - 

^ Re R Q 



F^gQg being the diagonal el- 



Here Ws s 

Rl e Rl <' 



ements of the Fock matrix F . 

~ S '^Ri g Rt lR2<y Rt]''^Rt + ^ ^Ri g R2<y^ ^Ri g R2 o- 
Rt 

comprises the negative of the Hartree-Fock potential and 
also the off-diagonal elements of the Fock matrix. Note 
that F is block-diagonal and thus the Wannier orbitals 
are still separated into occupied and virtual Wannier or- 
bitals. 

Using the inverse generalized Wannier transfor- 
mation — as is done throughout this article — 
corresponds to the use of the full Fock matrix as zeroth- 
order Hamiltonian. In addition, a perturbative expansion 
of the off-diagonal elements in the Fock matrix is applied 
in succession of the inverse Wannier transformation. This 
conduct is equivalent to employing the partitioning of the 
Fock matrix introduced in Eq. (|31|) . 

The one-electron interaction is mediated by the off- 
diagonal Fock matrix elements, and the two-electron in- 
teraction is represented by the two-electron integrals; let 
the one-electron interaction be treated using Feynman di- 
agrams up to m-th order and account for diagrams which 
exclusively involve two-electron interactions of highest 
order n. This separation facilitates to take higher-order 
terms into account for the slowlier convergent series; the 
resulting scheme is denoted by C0-ADC(r7i,n). 

The CO-ADC(2,2) approximation for Wannier orbitals 
is obtained similarly to the case of the C0-ADC(2) ap- 
proximation for Bloch orbitals H18|l and reads 



ui 

-R e ; <?i K 92 Q 93 
K+ 

gi K g2 a \ 
g{ «' 33' a' S3' /3' 

r+ 



^RggiK [32 a 93 /3] " a ^ , 

(5 91 K 92 o, 93 fi ; (Sgj « ~l" ^^93 
9/ So a' g^ /3' 



£91 ft) '^91 n '^92 a n-gs : 



= 



(32a) 
(32b) 

(32c) 
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The equations for U , K , and C are formally iden- 
tical to Eq. H^^2|l upon modifying the labels according 
to the adaptation of Eq. to M~g{k,uj) and replac- 
ing the occupation numbers in Eqs. H32a|l and l)32b|l 
by 'rtg'i a n-ij^ K n-gsX- Up to second order, the perturbative 
treatment of the off-diagonal elements of the Fock matrix 
does not introduce any new diagrams for the dynamic 
self-energy and thus K + C remains diagonal i^SiSi 
However, in a formulation which is invariant under uni- 
tary transformations in the occupied and virtual space, 
this is no longer the case because the full Fock matrix 
needs to be chosen as the zeroth-order Hamiltonian in 
place of 7?WF theniiiS^ 



IV. STATIC SELF-ENERGY 

A. Crystal momentum representation 

The static self-energy has not been determined so far. 
It is represented by all those diagrams in the diagram- 
matic expansion of the self-energy where the external 
points of the diagrams correspond to equal times^ In 
a strict second-order treatment of electron correlations, 
using Bloch orbitals l(2Jl, the static self-energy is zero be- 
cause the first static self-energy diagrams arise the ear- 



liest in third order. Moreover, it turns out that the 
diagrammatic series for the static self-energy does not 
converge reasonably in many casesiSi However, a self- 
consistent solution is possible^i*^ utilizingS2*SL2& 



k's] 



1 
2m 



§ Gsrik', u>) duj 



(33) 



here given in terms of Bloch orbitals. The contour in- 
tegration in Eq. H33|) runs along the real axis and closes 
in the upper complex w-plane, hence enclosing only the 
poles of the advanced Green's function G^{k',uj) in 
Eq. |6B)l . 

The self-consistent solution of Eq. (|33|l is computation- 
ally expensive. Yet a stable and efficient Dyson expansion 
method (DEM) for determining the static self-energy is 
devised in Ref. There, the first two terms of the 
Dyson expansion (|9bp are inserted into Eq. (|33|l . The 
term Wpq{k) in Eq. cancels the result of the con- 
tour integration over the free Green's function G^^{k', uj). 
Carrying out the contour integration over the product of 
the two free Green's functions and the static self-energy 
yields 



k p k q. k' s k' r k p k' r [k q k' s] ^ k' r k' s ' 

bpq{k) = y^,y^,Vkpk'r[kq k' .^] Qrsjk') , 



k' r k' s k' r k' 



s{k') - ^ ^ G°,(fc',c^)A/,,(fc',c^)G°,(fc',c^)dw 



(34a) 

(34b) 
(34c) 

(34d) 



The Qrs{k') can now be determined by inserting the spectral representation of the dynamic self-energy l|13(l into 
Eq. Ij34d|) . Carrying out the contour integration yields 



Qrs{k') = Q+{k') + Q;sik') , 



J2 m+'(")(fc')m+'(")*(fc') 

ne{N+l} 

nr, fir, 

k' r k' s 



k r k 



(er,~uj+{k')){ep-uj+(k')) 



k' r k' s 



Qrsik') 



r - ^k' s)i^k' r - ^t{k')) (Efc, r " ^k' s)i^k' s ~ ^rt (fc')) 

fir, fir, 

k' r k' s 



(35a) 



(35b) 



E m7^(")(fc')m;^(")*(fc') 
ie{N-i} 

nr, nr, 

k' r k' s 



{ej:,^-wn{k')){e^, ^-ujn (fc')) 
nr, nr, 

k' r k' s 



^ H _ . 

I 



(35c) 



In order to evaluate Eq. H35|) . from now on, we resort to the ADC form of the dynamic self-energy in terms 
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of Bloch orbitals (|17(l . By affixing a bar accent to the 
quantities U^ik), K^ik), and C^(k) in all following 
equations of this subsection, they become formally iden- 
tical to the ones obtained with the fully translational 
symmetry-adapted form H29() . 

The direct usage of Eq. I|35|) is in practice very time 
consuming because the K^{k') + C^{k') matrix has 
to be diagonalized fully H19() to obtain the pole posi- 
tions uinik') of the spectral representation ifT^ . For 
the C0-ADC(2) approximation of the dynamic self- 
energy in crystal momentum representation p8|l . how- 
ever, K^{k') + C^{k') remains diagonal. This leads 
to m^^"^*(fc') = Ut and w±(fc') = K^Ak') where 

k' r.n " 

I 



We have reduced the problem of determining 
the Qfg{k') to the problem of determining the V^{k'). 
Inspecting Eqs. H36|) and H38I) . we note that we have to 
solve for the Vj^{k') in the large 2pl/i-configuration space 
but only for the usually small set of all occupied orbitals. 
Conversely, the V^{k') have to be calculated in the small 
2/ilp-configuration space but for the usually large set of 
all virtual orbitals. 

The diagonal parts Sj^,^! — K^(k') in Eq. (|36|l are, 
at least, of a magnitude of twice the Hartree-Fock band 
gap, implying usually a diagonal dominance of the full 
matrices in Eq. H3t)|) . Therefore, a solution by Jacobi 
iterations'"'^'* is suggested: 

y±'(«)(^') = {eJ:,^-K^{k'))-^Ut{k'). 
Vr^'^"\k') = Vr^'^"\k') + {ef:,^l-K^{k'))-' (39) 
X C^{k')Vr^''-"-^\k') , 

which turns out to converge rapidly. 

The inhomogeneous linear system of equations H34afl 
can now be solved for I]^(A:) by a matrix inversion^: 

= (1 - A)-i5, (40) 



the n-th eigenvector of K^{k') + C^{k') is unity for a 
single 2plh- or 2/ilp-configuration, respectively, and zero 
otherwise. In general, the Qrs{k') are determined ef- 
ficiently from Eq. (|35(l by the inversion method or by a 
single- vector Lanczos diagonalization to circumvent a full 
diagonahzation of K^{k') + C^{k')^ 



The inversion method has been found to be more efh- 
cient£-° than the single- vector Lanczos diagonalization, so 
we concentrate on the former. We define auxiliary vec- 
tors {k') as the solution of the inhomogeneous systems 
of linear equations 



(36) 



(37) 



(38a) 
(38b) 



where hi, and Aij are composed by numerating the 
compound indices (p, q, k) and (s, r, k') in I]^(A:), b^{k), 
and Aj:^ k' sk'r ^"7 integer numbers i and j, respec- 
tively. Inspecting Eq. H34bp . we note that A couples 
only to the components of E^(fc') where the s,r indices 
denote a combination of particle-hole- or hole-particle- 
orbitals. As all quantities in Eq. (|40|l are Hermitian, we 
can further restrict the band indices to p < q. As soon 
as the number of fc-points becomes large, the solution of 
Eq. H40I) gets cumbersome as the dimension of the system 
of linear equations is given by the number of entries of 
the upper triangle (including the diagonal entries) of the 
static self-energy matrix times the number of fc-points. 
Thus an iterative linear equations solver has to be em- 
ployedi24 



B. Wannier representation 

The static self-energy can also be determined directly 
in Wannier representation. Having expressed the static 
self-energy in terms of Wannier orbitals, we obtain the 



(££,^1 - K±(fc') - C±(fc')) Vr^{k') = UHk') . 

Inserting the eigenvector matrix Y^{k') from Eq. H19|l into Eq. (|36|l . one obtains, by using formula (|20|l . 

y±(fc') = Y^ik') (££,^1 ~ fl±(fc'))-i™±t(fc') , 

which reveals the usefulness of the new vectors and allows us to rewrite Eq. (|35|l as follows: 

g+(fc') = -V+^ik')V+(k')n^^,^nj:,^ + {ef:,^-spJ-' 

X [c7+t(fc') V+ik') ng, ^ n-,, ^ - V+^k') U+ik') ng, ^ ng, J , 

g-,(fc') = V-\k')V-(k')n~^,^n^,^-iep^^Sf:,X' 

X [U-^ (fc') V- (fc') n^^, ^ ^ - V-Hk') LI~ (fc') n^, ^ nj:, J . 
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static self-energy in crystal momentum representation by 



5]e"=^E-(i?). (41) 

R 



To find an approximation for S ]^(i? ), we insert the in- 
verse Wannier transformation (|2rbp and its Hermitian 
conjugate into Eq. H34|l and arrive at 



i?' 9 
+ bga{R) , 



E-(i?') (42a) 



Q,„(i?') . (42b) 



Note 



that 



the 



translational 



symmetry 



of As 5 ^ ^ , s, could only be utilized once to 
remove one lattice summation such that two lattice 
summations (one over R' and one over g) show up in 
Eq. I|42|) while only one summation over k' is necessary 
in Eq. (|34|) . There is no first-order contribution to 
the static self-energy in Eq. H42() because the static 
self-energy used here evolves from an inverse Wannier 
transformation of the self-energy obtained in M0ller- 
Plesset partition£L22iS and a subsequent perturbative 
expansion with respect to the off-diagonal elements of 
the Fock matrix. 

For the remaining equations (|34bp and (|34dl) . we as- 
sume the partition (|31|l of the Hamiltonian in Wannier 
representation, to obtain 



A- - - - - 

g R(T. g V g+R' r 



IrviR ) — Qj^r ^ g 



^Og g+R' T [Ra gv] 

^ i^g+R'r~^av)'' (43a) 

1 



27ri 



Gl,{G,Lo) (43b) 



with the free Green's function in Wannier representation 
given by 



GrviR' ^ R^ 



^Rt. R' V 



(44) 

To determine the Qtv{R'), we introduce the spectral 
representation of the dynamic self-energy in Wannier rep- 
resentation [the analog of Eq. (|13|l ] 



M,r{R\oj) = 



E 

ne{N+l} 



mj'(")(0) m;^^"^*(i?') 



+,(«)*/ 



UJ ■ 



IT] 



ie{N-i} 



(0)m7'^"'*(i?') 



IT] 



(45) 



M+,{R',iu) + M-,{R\u;) 



which is obtained from the ADC form of the dynamic self- 
energy in Wannier representation (|23|l by letting R ^ 

and diagonalizing K + C analogously to Eq. H19|l . 
Since translational symmetry cannot be exploited in 
Eq. as R' — in contrast to k — is not a good quantum 
number, this leads to the redundancies that have already 
been mentioned in conjunction with the superccU form of 
the dynamic self-energy in Eq. (|25|) ; namely, the full con- 
figuration space of 2plh- and 2/ilp-configurations from 
the A^o unit cells of the crystal is necessary here. How- 
ever, the cutoff criterion of Sec. I VII nevertheless, can be 
applied. 

The Qtv{R') can now be determined by inserting 
Eq. (|45|) into Eq. Ij43b|l . Carrying out the contour in- 
tegration yields 



Qr,{R') - Qt,.{R') + Qr.{R') , 

QUR') ^ rhP^"HO)rh+'^"^*{R') 

ne{N+l} 



-'^R'r'^Ov 



r - ^Ou)i^R' r " ) R' r ^ '^0 ^,)i£o v ~ ) 

ie{N-i} 



(0)m-^(")*(i?') 



''R' T "-Q u 



i^R' r - ^dJi^R' r " ) R' r " ^Ov)i^Ov " ) 

We proceed as before by defining the auxiliary vectors 



(46a) 



(46b) 



(46c) 



(47) 
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setting Uf{R) = Uj^ , and rewrite Eq. (gSJ: 



)UR') = -y+t(0)y+(i?')"^,,no, + (£^,,-£oJ-^ 

X [i7+t (0) V+ {R') n^, , nj^ - l7+t (0) «7+ (i?') n^, , J , (48a) 

- K-t(O) . ^0. - i^R' r - ^oX' 

X [c7-t (0) V- {R') n^, ^ ng, - V'^ (0) C/- (i?') n^, ^ fig J . (48b) 



The algorithm for Jacobi iterations now reads 

V,^'^'\r') = iej^,l-K^)-'U^iR'), 
l/±'(")(i?') = F±^(°)(i?') + (e^,,l--^*)-^ (49) 
xC^l?±^^"-^)(i?') , 

and the inhomogcneous hnear system of equations H42a|l 
can be solved for E^(i?) by a matrix inversionSi: 

t ^{i-A)-'^b, (50) 

where a check accent is affixed to indicate that Wannier 
orbitals rather than Bloch orbitals are used. We employ a 
mapping of the compound indices {g, a, R) and (v, r, R') 
to the integer numbers i and j, respectively, for the quan- 
tities Ef, hi, and Aij. Note that S°°(^) is not Hcrmi- 
tian. Hence, the full static self-energy matrix has to be 

included in E explicitly for several lattice vectors R. 
However, the property t(^) = S°°(-^) holds which 
can be utilized to reduce the number of equations in for- 
mula (|5U|) . Nevertheless, in many cases, the linear system 
of equations needs to be solved iterativelyiSi 



self-consistently, i.e., such that the energy u entering the 
dynamic self-energy M{k,oj) is identical to the resulting 
eigenvalue lu. 

To avoid the unitary matrices V{{k) and W^(fc) in 
Eqs. (|29a|l and H41() arising in the following expressions, 
we multiply Eq. (|51|l from the left with lA^k) and in- 
sert ti{k)ti'^ (k) — 1 between the right parenthesis and 
the eigenvector xq on the left-hand side of the formula. 
Defining the translational symmetry-adapted Fock ma- 
trix Fga{k) — ^ e'*^^ Fgcr(^), one arrives at 
R 

iF{k) + -E"°{k) + M{k,uj))xGik,Lu) ^ LuxGik,Lu) , (52) 

where the new quantities Fga{k), E^(fc), Mga-{k,uj), 
and XQ^g{k, uj) are related to the old ones epq{k), E^(fc), 
A/pq(fc,cj), and XG,p{k,uj) by e{k) = U{k) F{k)U\k), 
S°°(fc) = U(k)i:°°(k)U''{k), M(fc,w) 
U{%)M{k,uj)U\k), and ±G{k,uj) ^U{k) XG{k,uj). 

Inserting the ADC form of the dynamic self- 
energy H29a|) into Eq. H52|l and defining 



V. DYSON EQUATION 

Having determined approximations for the static and 
the dynamic self-energy in terms of Wannier orbitals in 
Eqs. and (j^ . respectively, we can finally determine 
the positions and strengths of the poles of the one-particle 
Green's function from Eq. (|10l) . [The pole positions and 
strengths of the Green's function, using Bloch orbitals to 
represent the self-energy (I34f) and H17|) . can be obtained 
following a nearly identical line of argument.] To this 
end, the small Hermitian eigenvalue problem 

{e{k) + Y.°^{k) + M{k,Lo))±G{k,Lo)=Lo±G{k,to) , (51) 

with e{k) being the diagonal matrix of Bloch orbital ener- 
gies and XG{k, uj) denoting eigenvectors, has to be solved 



f ± (fc, uj) = {uj1-K^ {k) ^C^{k))-^ (fc) XG {k, oj) 

(53) 

yields the following form of the eigenvalue problem l|52|l : 

{F{k) + i:'^{k))xG{k, u;) + U+\k) x+{k, u;) ^^^^ 
+ U \k) x^ {k,uj) — UJ XG{k,uj) . 

Rewriting Eq. I|53(l as 

U^ik)xG{k,uj) + {K^{k) + C^$))x^{k,Lu) 
— ujx^(k,uj) , 

it becomes evident that Eqs. (|54|l and H55|l can be recast 
as a joint Hermitian eigenvalue problem 



B{k)X{k) = X{k)E{k), X\k) X{k) = 1 , 
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B{k) 



f F{k) + tr{k) 

\ u-{k) 



U^\k) 
K+{k) + C+{k) 



u-\k) \ 



K'{k) + C'{k) ] 



(56) 



The new matrix B(k) is called the hand structure matrix. 
It has to be diagonalized for several fc-points yielding 
eigenvalues e„(fc) = {E{k))nn and eigenvectors X„(fc) = 
(fG(fc,e„(fc))T, a;+(fc,e„(fc))T, e„(fc))T)T. For- 

mula (|56|l is similar to the result for moleculesi^ 
One can envisage that the diagonalization of the sub- 
blocks K^{k) + C^{k) in B{k) is carried out first by 
solving Eq. (|19|l : afterwards, the full eigenvalue prob- 
lem (|56ll is treated with all subblocks modified accord- 
ingly.^'^'*'* Therefore, via Eq. H15|) . infinitely many proper 
self-energy diagrams are summed before the result is 
put in the Dyson equation (j^ to sum all improper self- 
energy diagrams that derive from those contained in the 
ADC form. 

The spectral representation of the one-particle Green's 
function © reads, in terms of the eigenpairs of the band 
structure matrix (|56|l . 




(57) 

where x^p^\k) = {X{k))pn denotes the transition ampli- 
tude of the n-th state. It is given by either the x'p'\k) 
or the j/p"-* (k) in Eq. whether it belongs to an — 1 
or an -I- 1 particle state. 

In some cases, the external Wannier orbitals of the 
dynamic self-energy are restricted to a single cell of the 
crystal, i.e., M{k,uj) w M{R = 0, w), the so-called one- 
lattice-site approximation for the dynamic self-energy. 
Note, however, that the 2plh- and 2/ilp-configurations 
in neighboring unit cells are not entirely excluded. The 
fc-dependence, in this approximation, is solely due to the 
one-particle matrix F(fc) + S°°(fc). The static self-energy 
can also be evaluated in one-lattice-site approximation. 
Then the fc-dependence is exclusively due to F{k) and 
this approximation to the CO- ADC equations becomes 
similar to the dynamical mean-field theory^ 

VI. CONFIGURATION SELECTION 

All equations derived so far contain certain infinite lat- 
tice sums that run over the entire crystal forming the con- 
figuration space. To obtain a meaningful approximation 
of the one-particle Green's function, the convergence of 
those lattice sums must be granted. Sun and BartletliS^ 
study the convergence properties of the lattice sums in 
several quantum chemical correlation methods for ground 
states and exited states of crystals, based on the M0ller- 



Plesset partition of the Hamiltonian Following the 
reasoning of Ref. 95, the lattice sums in CO- ADC in crys- 
tal momentum representation converge. Moreover, the 
equations in local representation also converge as they 
evolve from the former by means of an inverse Wannier 
transformation (|21b|l . 

Although the convergence of lattice sums is granted 
principally, we have to devise an algorithm for their 
proper truncationj^ i.e., a dynamical building of the con- 
figuration space is required which means to meet a cho- 
sen accuracy in the lattice sums for a given crystal. The 
appropriate truncation of lattice sums is an essential in- 
gredient of all ab initio methods for crystals and corre- 
sponds to a configuration selection procedure which en- 
sures that the configuration space is sufficiently large for 
the desired accuracy of the band structure but still suf- 
ficiently small to be tractable on present-day computers. 
In molecular physics, configuration selection has been in- 
troduced in the context of the configuration interaction 
method— which suffers from an exponentially growing 
configuration space with respect to the number of atoms 
in molecules. The techniques discussed in Ref. 96' are not 
applicable here, as they are designed with the configura- 
tion interaction method in mind, but there are analogs 
for crystals which are similar in spirit. 

The incremental scheme for ground states of crys- 
tala'*i^i^^i'^^i^*iQ^i^°°i^°^ is based on Wannier orbitals. The 
total configuration space of a crystal is partitioned into 
certain subsets which are used in a subsequent calcula- 
tion of the correlation energy to yield the so-called energy 
increments. Frequently, explicit manual configuration se- 
lection is applied. This allows the partition of the config- 
uration space to be chosen following chemical intuition. 
Yet the situation is much more cumbersome for band 
structures because the number of distinct excited states 
grows with the number of correlated electrons, i.e., the 
number of Bloch orbitals and Wannier orbitals which are 
considered in the sums in Eqs. (|17|l and H29(l . respectively. 
In the local Hamiltonian approach, the number of states 
to be calculated is restricted to a fixed number which 
is treated in an incremental ■\vayi^2i^2i24i2Sii£ Albrecht et 
alm^i^ ensure the convergence of an incremental series 
for the self-energy for w = 0, a criterion which is inde- 
pendent of the number of states actually described. 

We devise here a configuration selection procedure 
that is perfectly adapted to the structure of the CO- 
ADC scheme. For the equations in terms of Bloch or- 
bitals (|34|l and (|17|) , fc-points are chosen harnessing Born 
von Karman boundary conditions which lead to a net of 
equidistant fc-points."'^'^ In this case, configuration selec- 
tion means choosing a sufficiently dense net in conjunc- 
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tion with a cutoff criterion which is similar to the one for 
CO- ADC in terms of Wannier orbitals (gTJ and (|^ dis- 
cussed in the next paragraph. Yet beforehand, one should 
interpret the impact of the number of fc-points used for 
Brillouin zone integration. Using a given net of A^o k- 
points to carry out a Wannier transformation 121a() , CO- 
ADC in terms of Bloch orbitals can be immediately rep- 
resented and analyzed in terms of the supercell form (|25|l . 
However, one should emphasize that the redundancies of 
this representation — which is utilized here for interpre- 
tation only — are not present in the Bloch orbital formu- 
lation of CO- ADC. The supercell form affords a tight- 
binding interpretation: within the supercell, the number 
of distinct neighbors to the origin cell is restricted by 
the volume of the supercell. Distinct interaction terms 
among the cells, i.e., the Fock matrix and two-electron in- 
tegrals, are thus also resticted and the number of fc-points 
employed can thus be understood to imply the number 
of nearest-neighbor cells treated in the interaction terms. 
These arguments gene ralize the tight-binding arguments 
given, e.g., in Ref. 'SCf, in conjunction with a single Bril- 
louin zone integration of a crystal-momentum-dependent 
function to the multi-dimensional case. 

Configuration selection in crystal momentum represen- 
tation is not intuitive and does not allow a fine-grained 
selection of configurations within unit cells and between 
an origin cell and its neighbor cells. For crystals with 
large unit cells this is a significant restriction. To obtain 
a suitable cutoff criterion for CO- ADC in terms of Wan- 
nier orbitals H41I) and (|29|l . we evaluate the second-order 
diagram of the retarded dynamic self-energjiM*^2iS and 
examine the summand therein: 
l/- V* 

e 31 f= [S2 a 33 /3] Ra giK [pa a gs P] 

n-gi K ^T-ga a ng^ p . 

LU fcgaa ^g^P'^giK, 

(58) 

The most delocalized occupied orbitals are the valence 
orbitals. Therefore, the occupied valence orbitals cou- 
ple to the most important 2pl/i-configurations. Hence, 
uj can be assumed to be of the order of the band 
gap. The denominator in Eq. (|58|1 is regarded to be 
constant to a reasonable degree. It is thus consid- 
ered as being part of the cutoff threshold. Therefore, 

if ^0(?Sl«;[92og3/3]'^j^agiK[g2aS3/3] '^^^ ""S^ « ""93 IS 

above a certain cutoff threshold for some combination 
oi Q Q and Ra, we include the corresponding 2pl/i-test- 
configuration [gi k, 32 a, 53 0) in the configuration space. 
If the product is small for all external Wannier orbitals 
of the self-energy, {gi k, g2 a, gz (i) is neglected completely 
by setting all two-electron integrals containing this con- 
figuration exactly to zero. By this, only a finite range 
of the residual Coulomb interaction, which is not treated 
in Hartree-Fock approximation, is explored. The same 
arguments hold for the selection of 2/ilp-configurations 
for the evaluation of the advanced dynamic self-energy, 
if Eq. (|58|l is adapted in the same way as Eq. (|32|l is 
adapted to C7 , K , and C . 

Convergence of the configuration space, i.e., the num- 



ber of unit cells contributing configurations to be consid- 
ered in the cutoff criterion (|58|l . is expected to be suf- 
ficiently quick. Van dcr Waals dispersion interactions 
which are present both in the ground state of the N par- 
ticle system and the states of the A^ ± 1 particle system 
are nearly equal and thus almost cancel in band structure 
calculations. However, the effect of the extra Coulomb 
charge that arises in A^ ± 1 particle states in contrast 
to the A^ particle ground state is long range. The poten- 
tial of an extra point charge falls off like V{r) = ±l/(e r), 
where e is the dielectric constant of the solid. This causes 
an overall polarization of the crystal which can be ac- 
counted for by adjusting the chemical potential /i^ of 
the added or removed electron accordingly by using a 
continuum approximation4i^2i^2i^ This procedure is not 
specific to the CO- ADC theory advocated here to treat 
the short-range interactions, but is required generally. 

Having introduced a convenient method to truncate 
the configuration space for CO- ADC, we have to elu- 
cidate the physical impact of configuration selection on 
band structures. To this end, we consider the fully 
translational symmetry adapted form (|29(l . The mod- 
ified coupling amplitudes l|29b() in the band structure 
matrix H56|l of a crystal carry two independent crystal 
lattice vectors g{ and If 9i ^^'^ 92 sufficiently 
far away from the origin cell, all two-electron matrix el- 
ements (j32all which contribute to the modified coupling 
amplitude [7^ , , „,(k) fall for all a below the con- 

figuration selection threshold and are therefore zero. Be- 
cause the modified interaction matrices (|32c|l vanish in 
the CO-ADC(2,2) approximation and Eq. H29c|l thus de- 
notes a diagonal matrix, the eigenvalue associated with 
the configuration {g[ g^ k! o! /?') is simply given by the di- 
agonal element {K {k))gi ^, p, . g> g> ^, p, , with an 
eigenvector being unity on position gl g!^ n' a' j3' and zero 
elsewhere. Consequently, such an eigenstate of B{k) is 
removed by deleting the column and the row g[ 52' /?' 
in -B(A:) because the transition amplitude lO vanishes 
and thus it does not contribute to the spectral represen- 
tration of the one-particle Green's function 1)57(1 . 

A further consequence of configuration selection in- 
volves the lattice summations over R in Eqs. I|29b|l and 
(I29c(l . In the first place, we consider a crystal with a 
macroscopic lattice constant which consists of Aq unit 
cells. As a result, only Fock and two-electron matrix 
elements that involve Wannier orbitals from the origin 
cell are selected. Therefore, in Eqs. (|29b|l and H29c() . the 
lattice sums run over only the single term for i? 
with gi = g2 — and the resulting band structure ma- 
trix (|56|l is consequently independent of fc. As the Bril- 
louin zone contains A^o crystal momenta, we obtain a 
A^Q-fold degenerate spectrum similarly to simple tight- 
binding modelsi In typical crystals, interactions with 
neighboring unit cells are important. Yet a sufficiently 
large supercell which consist of no unit cells can be chosen 
such that interactions between supercells are negligible 
leading to a Nq/uq fold degenerate spectrum. In other 
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words, a cutoff threshold, to be used for the selection of 
Fock matrix elements and 2plh- and 2/ilp-configurations, 
implies a definition of degeneracy of the physical states 
in our model of the crystal. 

In order to selectively diagonalize the band structure 
matrix, one uses an iterative eigenvalue solver, e.g., a 
block-Lanczos algorithm^^^iSiiSSiiS^ that is capable of 
exploiting the sparsity of B(k). The most expensive 
step of the block-Lanczos algorithm is a matrix times 
vector producljiii between B{k) and Lanczos vectors 
which determines the overall performance of the eigen- 
value solver iSii2iii2ii24 To investigate the scaling of a se- 
lective computation of eigenpairs, we consider a supercell 
consisting of two unit cells of an original lattice. Let us 
assume a crystal with a macroscopic lattice constant to 
investigate the asymptotic scaling behavior of the prob- 
lem. Then, the configuration selection method of the 
previous paragraphs selects only configurations local to 
the individual unit cells. Hence, the resulting configura- 
tion space scales linearly with the system size, i.e., the 
number of atoms per unit cell. This is a necessary con- 
dition for the selective diagonalization of the band struc- 
ture matrix (|56|l to scale linearly as well. The total Fock 
matrix of the supercell also is block-diagonal. Similarly, 
the modified coupling amplitudes (|29b|) only contain two 
nonzero blocks which describe the coupling of 2plh- and 
2/ilp-configurations to the Fock matrix of each of the two 
constituting unit cells of the supercell such that the total 
band structure matrix decomposes into two subproblems 
that can be solved independently. A single matrix times 
vector product for the supercell requires twice as many 
floating point operations than are needed for a matrix 
times vector product for one of the two unit cells in the 
supercell. Hence, the computation of matrix times vec- 
tor products scales linearly with the system sizCi^i^ Yet 
doubling the system size also usually means that we are 
interested in twice as many excited states. Hence, the 
overall effort to determine all excited states of a crys- 
tal in a given energy range scales quadratically. It is a 
quadratically scaling problem, where linear scaling can 
only be achieved by an a priori restriction to a few ex- 
cited states of the system. 

For typical crystals, the Fock matrix of a supercell is 
not block-diagonal which implies another doubling of the 
number of floating point operations upon doubling the 
system size. Yet this factor cancels as the number of 
fc-points needed for a given accuracy of the integration 
over the Brillouin zone of the supercell is halved because 
the volume of its Brillouin zone is half the volume of the 
Brillouin zone which corresponds to the original crystal 
latticed 



VII. LITHIUM FLUORIDE CRYSTAL 

Lithium fluoride crystallizes in a face-centered- 
cubic (fee) rocksalt structure described by the space 
group FmZm as shown in Fig. ^ The crystal lattice 




FIG. 1: (Color) Structure of a LiF crystal. The lithium and 
fluorine atoms are represented by solid yellow and hatched 
green spheres, respectively. 



has a two-atomic basis Li(0,0,0) and F(i, i, i), given in 
units of the lattice constan'feiiS a = 3.990 A. To inves- 
tigate a LiF crystal, we have chosen the simplest pos- 
sible ah initio description. The main purpose here is 
to demonstrate the feasibility of the CO- ADC formulas 
for a three-dimensional crystalline solid. Thus, only a 
minimal Gaussian ST0-6G basis se1ii2^*i2& is used which 
describes the lithium atom by a (12s) / [2s] contraction 
and the fluorine atom by (12s 6p) / [2s lp\. Each shell of 
the ST0-6G basis set is constructed by fltting six Gaus- 
sian functions to each occupied orbital of the isolated 
atoms. The 2s shell of a lithium atom is rather diffuse 
due to the single outer valence electron. As the lithium 
atoms are ionized in LiF crystals, we remove the two 
most diffuse Gaussian functions in the Li 2s contraction 
because they do not reflect the physical situation of a 
compact Li+ ion, and arrive at a (10s) / [2s] contraction 
scheme. 

The Hartree-Fock equations are solved self-consistently 
by means of the wannier program^Sii^ which directly 
yields Wannier orbitals. A finite cluster of unit cells is 
utilized as support for the Wannier orbitals of the origin 
cell which, in our case, consist of up to third-nearest- 
neighbor cells (43 unit cells alltogether) . As the pro- 
gram determines only occupied Wannier orbitals, crystal 
projected atomic orbitals (crystal PAOs)^i4i2iiiS£ have 
been devised on the basis of the projected atomic orbitals 
introduced by Saeb0 and Pulajii22ii2S and also used by 
Hampel and Werneriiifl The resulting virtual functions 
are subsequently Lowdin orthonormalized^Siii'2L2& For 
the CO-ADC(2,2) method and the Fock 

matrix and the two-electron integrals in Wannier rep- 
resentation, entering the Hamiltonian (|31|l . need to be 
obtained by a transformation of the corresponding quan- 
tities in terms of Gaussian basis functions. This trans- 
formation is presently the bottleneck in practical com- 
putations since the wannier prograrr42ii^ carries it out 
for clusters only and neither uses translational nor point 
group symmetry. At present, this permits only the use of 
a minimal basis set for a three-dimensional LiF crystal. 
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16.24 


16.24 


4.95 
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0.41 


16.16 


15.74 


4.80 


13 


2.43 


15.69 


13.26 


4.78 


19 


2.89 


15.36 


12.47 


4.75 



TABLE I: Convergence of the top of the valence bands i5top,v, 
the bottom of the conduction band i?bottom,c, the fundamen- 
tal band gap -Egap, and the bandwidth of the F2p valence 
band complex AEf2p of a LiF crystal. "Cells" designates 
the number of unit cells taken into account in CO-ADC(2,2) 
calculations of the quasiparticle band structure where zero 
refers to the Hartree-Fock result. Unity denotes the inclu- 
sion of the 2plh- and 2/!,lp-configurations from the origin cell 
only in correlation calculations. 13 and 19 indicate the ad- 
ditional inclusion of configurations from nearest- and next- 
nearest-neighbor cells, respectively. For a LiF crystal, Etop.v 
and -Bbottom.c are both situated at the F point. All data are 
given in electronvolts. 

The subsequent correlation calculations with the CO- ADC 
program^iiiii are very efficient and as rapid as typical 
molecular ADC calculations. 

The Hartree-Fock band structure that results from this 
minimal basis is shown in Fig. |21 All band energies are 
given with respect to the top of the valence bands which 
is located at the T point and is set to zero. Below zero 
there is a complex of three F 2p valence bands shown as a 
closeup in Fig. O They originate from the three F 2p va- 
lence energy levels of the LiF molecule around 12 eV 
where the lowest ionization potential of the molecule is 
doubly degenerate^i Upon crystallization, this degener- 
acy is lifted due to interactions with neighboring atoms 
and, except for the high symmetry lines F-L and F-X, 
one observes three distinct F 2p valence bands. As a mini- 
mal basis set is utilized, there is only a single conduction 
band of mainly Li 2s character. A slight admixture of 
F 2p character to this band is observed in the W~K panel 
by comparing with the band structure resulting from a 
larger basis set, e.g., the one of RcLIsS where a conduc- 
tion band of F 2p character mixes with the Li 2s band. 
Nevertheless, the single conduction band represents the 
energetically lower edge of the conduction-band complex 
well. LiF has a direct band gap; i.e., the maximum of the 
F 2p valence bands and the minimum of the Li 2 s conduc- 
tion band are located at the same crystal momentum, the 
F point, here. 

In order to investigate electron correlations, we start 
from a single unit cell, to form the 2plh- and 2hlp- 
configurations entering Eq. (|29|l for M^{k,uj), respec- 
tively. Successively, configurations which contain Wan- 
nier orbitals in nearest- and next-nearest-ncighbor cells 
are considered. This leads to three distinct quasiparticle 
band structures whose convergence with respect to the 
number of unit cells where electron correlations are re- 
garded explicitly is monitored for the following key quan- 
tities: the top of the valence bands, the bottom of the 
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FIG. 2: (Color) Band structure of a LiF crystal. Hartree- 
Fock bands are given by solid green lines. The dashed red 
lines depict the CO-ADC(2,2) quasiparticle bands which are 
determined by accounting for 2plh- and 2/ilp-configurations 
involving Wannier orbitals in the origin cell and the nearest- 
and next-nearest-neighbor cells. 



conduction band, and the band gap and the width of the 
F 2p valence-band complex [Tab.|l]. The results are com- 
pared with the plain Hartree-Fock data (with configura- 
tions out of "0" unit cells). Taking into account electron 
correlations of a single unit cell causes a slight upwards 
shift of the valence bands by 0.41 eV while the conduc- 
tion band essentially remains unchanged. Upon inclusion 
of configurations in the nearest-neighbor cells, the band 
structure changes drastically, leading to a significant re- 
duction of the band gap by 2.48 eV. Accounting for con- 
figurations which extend to next-nearest-neighbor cells 
has three times less impact on the band gap. Therewith, 
it can clearly be seen that the major contributions of elec- 
tron correlations have been covered and that the effect 
of configurations which extend to Wannier orbitals be- 
yond next-nearest-neighbor cells will yield much less sig- 
nificant changes of the quantities than those discernible 
in Tab. ^ As pointed out in Sec. IVII the effect of the 
more distant unit cells can be determined using a con- 
tinuum approximation. ^■^^•■^^■^^ The observed correlation 
corrections in Tab. ^ with a growing number of neigh- 
boring unit cells result primarily from the formation of 
a quasiparticle and the polarization of the surrounding 
by the Coulomb charge of the added particle or holei^ 
Note that, additionally, the beneficial basis set exten- 
sion (BSE )-t5fU 2 effect is noticeable in solids which leads 
to an improved description of electron correlations with 
an increasing number of unit cells taken into account. 

The quasiparticle band structure of LiF is displayed 
in Fig. 121 for the F 2p valence bands and the Li 2s con- 
duction band. The former bands are additionally dis- 
played on an enlarged scale in Fig. 13 The data are 
taken from the computation with configurations from all- 
together 19 unit cells. Valence and conduction bands do 
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FIG. 3: (Color) Valence band structure of a LiF crystal. Zoom 
into Fig. 121 



not shift by the same amount, upwards and downwards, 
respectively. Instead, the former bands arc much more 
strongly influenced by electron correlations than the lat- 
ter bands. This is due to the fact that the occupied 
Wannier orbitals are spatially less extended than the vir- 
tual Wannier orbitals. For electrons in compact orbitals, 
the Hartree-Fock approximation becomes progressively 
worse due to an insufficient description of the Coulomb 
holc.^ Although the quasiparticle bands are considerably 
shifted with respect to the corresponding Hartree-Fock 
bands, they essentially keep their width. We also observe 
this effect for (IIF)oo chains. This is in contrast to the 
significant reduction of bandwidths observed for quasi- 
particle band structures of covalently bonded polymers 
like trans-polyacetylena^ and covalently bonded crystals 
like diamond, '^^•■^•^^^ silicon^Si^, and germanium. '^'^ 

The ab initio description of a LiF crystal using a min- 
imal basis set provides valuable insights. However, we 
cannot expect full quantitative agreement neither at the 
Hartree-Fock nor at the correlation level. The accuracy 
of our approach can be determined by comparing with 
experimental and theoretical data from related studies. 
The fundamental band gap is accessible by photoelec- 
tron spectroscopy and has been measured to be 14.1 eV. 
Another experimentally accessible quantity, the width of 
the F 2p valence-band complex, was determined to lie 
in the range^S of 3.5-6 eV and thus is inconclusive, un- 
fortunately. With our ST0-6G-like basis set, we find 
a Hartree-Fock band gap of 16.24 eV [Tab.^ which de- 
viates considerably from the value of 22.7 eV communi- 
cated by Kuna£& or the value of 22.4 eV reported by Al- 
brecht^ who both use larger basis sets. For the width of 
the F 2p valence-band complex, we find 4.95 eV [Tab. 
in the Hartree-Fock approximation which is appreciably 
larger than the value of 3 eV from Kunz^*' and the value 
of 3.37 eV from Albrecht3 

Including electron correlations, the quasiparticle band 
gap reduces to 12.47 eV [Tab. ^ and thus falls short 



by 1.5 eV compared to the experimental value. Al- 
brechtS finds 13.5 eV, and calculations of Kunz^i yield a 
theoretical band gap of 14.0 eV, in good agreement with 
the recent measurements of 14.1 eV. For the ¥2p va- 
lence band complex, we obtain a width of 4.75 eV [Tab.|Il 
which is nearly unchanged compared to the Hartree-Fock 
value, a fact that has also been found in previous studies 
of bulk LiF.^^'^" Consequently, our value of the width re- 
mains much higher than the value of 3.1 eV from Kuna2£ 
or the value of 3.40 eV from Albrechti2£ 

The above results show that the CO- ADC theory can 
be applied without problems to compute the energy 
bands of three-dimensional crystalline solids. What re- 
mains to be improved considerably is the transformation 
of the Fock matrix and the two-electron integrals from 
the Gaussian basis set representation to a representation 
in terms of Wannier orbitals. 



VIII. CONCLUSION 

In this article, a general local orbital ah initio Green's 
function method is devised for band structures. The 
band structure is given by the pole positions of the 
one-particle Green's function in terms of Bloch or- 
bitals which in turn is evaluated efficiently utilizing 
Dyson's equation. To approximate the self-energy up 
to n-th order, we devise a crystal orbital formula- 
tion of the well-established algebraic diagrammatic con- 
struction (ADC)^^i^^i^ which is termed crystal orbital 
ADC (CO-ADC).Sa The pole search of the one-particle 
Green's function is recast into a Hermitian eigenvalue 
problem which is a numerically stable and efficient for- 
mulation, permitting us to explore strong correlations 
that occur for the energetically lower lying bands of cvys- 
tals.^'^^SiiLi2ii^ii4 Wannier orbitals are used in internal 
summations in the self-energy and allow us to exploit the 
fact that electron correlations are predominantly local. 

The lattice summations in the CO-ADC equations 
must be truncated to render the problem tractable. To 
this end, we devise a configuration selection procedure 
which can equally well be used in conjunction with the 
calculation of ionization potentials and electron affinities 
of large molecules to speed up computations. Configura- 
tion selection is shown to lead to a definition of degener- 
acy among the states of a crystal. 

The derivation of the local orbital CO-ADC theory 
sets out from the equations in terms of Bloch orbitals. 
Afterwards internal indices are transformed to a Wan- 
nier representation. We consider this line of argument to 
be compelling due to the close analogy of the equations 
in crystal momentum representation to the equations of 
molecular physics. Alternatively, the derivation can be 
conducted by starting from CO-ADC in terms of Wan- 
nier orbitals. A transformation of the equations to crystal 
momentum representation is then carried out by utilizing 
the Wannier transformation H21afl . 

The CO- ADC computer program^iiiii has been devel- 
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oped for computations based on the CO-ADC(2,2) ap- 
proximation. It requires the Fock matrix and the two- 
electron integrals in local representation which are ob- 
tained by an integral transformation using the Wannier 
orbitals from a Hartree-Fock calculation with the ab ini- 
tio program wANNlERi^Si^ At present, this transforma- 
tion poses a bottleneck because the current implementa- 
tion has been designed for molecular clusters and does 
not exploit any symmetries. Therefore, only a minimal 
basis set is used to study a lithium fluoride crystal in 
order to demonstrate that the method is working. A 
fundamental band gap of 12.47 eV and a bandwidth of 
the F2p valence-band complex of 4.75 eV is obtained. 
These values are to compare with the experimental data 
of 14.1 eV and 3.5-6 eV and the theoretical results of 13.5, 
14.0 eV, and 3.1, 3.4 eV for the band gap and the band- 
width, respectively. Increasing the size of the basis set 
will improve the CO-ADC(2,2) results. Yet it is demon- 
strated that the CO- ADC theory is well suited for tack- 
ling the problem of calculating energy bands of solids 
with controlled approximations. 



We would like to point out that the block of 2plh- 
configurations in the band structure matrix (|56|) is much 
larger than the block of 2/ilp-configurations. In order 
to reduce the computational effort, the advanced part 
and the retarded part of the one-particle Green's func- 
tion can be independently evaluated diagrammatically. 
This leads to a different, so-called non-Dyson CO- ADC 
scheme, which decouples the computation of occupied 
and virtual bands completeljiSiiH and results in two in- 
dependent Hermitian matrix eigenvalue problems to be 
solved. 
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